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We address the problem of interference management and power control in terms of maximization 
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of a general utility function. For the utility functions under consideration, we propose a power control 
algorithm based on a fixed-point iteration; further, we prove local convergence of the algorithm in 
the neighborhood of the optimal power vector. Our algorithm has several benefits over the previously 
studied works in the literature: first, the algorithm can be applied to problems other than network utility 
maximization (NUM), e.g., power control in a relay network; second, for a network with N wireless 
transmitters, the computational complexity of the proposed algorithm is 0(N 2 ) calculations per iteration 
(significantly smaller than the 0(N 3 ) calculations for Newton's iterations or gradient descent approaches). 
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f-*- , Furthermore, the algorithm converges very fast (usually in less than 15 iterations), and in particular, if 
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■<^j- . initialized close to the optimal solution, the convergence speed is much faster. This suggests the potential 



(N 
(N 



of tracking variations in slowly fading channels. Finally, when implemented in a distributed fashion, the 
algorithm attains the optimal power vector with a signaling/computational complexity of only O(N) at 
each node. 



^ " I. Introduction 



Interference control is one of the major problems of today's dense wireless networks; indeed, the 
latest releases of the Long Term Evolution (LTE) and 802.16 (WiMax) standards take measures to 
control the interference environment in dense network implementations (H, |2). Intra-cell interference, 
caused by other cellular users within the cell, occurs in cellular networks with non-orthogonal channel 
assignment within a single cell, e.g., code division multiple access (CDMA) networks. In more recent 
cellular network standards, this problem is resolved through orthogonal channel access assignments, 
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e.g., orthogonal frequency division multiple access (OFDMA) lU, |2). Inter-cell interference, i.e., the 
interference received from adjacent cells, is a significant problem even in modern cellular networks 
with intra-cell interference management. Essentially, the problem of designing practical and efficient 
interference control techniques has become extremely important in the design of next generation wireless 
cellular networks. 

A. Previous Work 

The power control problem to reduce co-channel interference was proposed as a linear programming 
problem as early as 1964 0, wherein the authors considered solving for the minimum transmission power 
given a required signal to noise and interference ratio (SINR). In the years that followed, researchers 
investigated power control in different applications. Most of the work up to the early 1990s focused on 
techniques to maintain the received power at a constant level (H (e.g., to tackle the near-far problem in 
CDMA networks). The work of Zander [15 ] and, soon after, the distributed power control (DPC) algorithm 
by Foschini et al. O were early solutions for SINR assignment in cellular networks. These preliminary 
works consider the problem of setting equal and arbitrary 01 target SINRs for wireless users in 
cellular networks. Although these solutions efficiently solve the power control problem in distributed 
cellular networks, they do not consider if the assumed initial SINR target is feasible. Moreover, the 
assumption of target SINR levels is mainly applicable to voice networks wherein a minimum SINR is 
required to maintain communications. In recent cellular network applications, more sophisticated power 
assignment techniques are necessary to ensure adequate quality of service (QoS) throughout the network. 

One body of research work following @ and |6l emphasize the distributed aspects of the power 
control algorithm at the cost of optimality of the solution (e.g., see Q-Q and the reference therein). 
For instance, the work in Q proposes a variant of DPC and proves its convergence. In their algorithm, 
the wireless nodes update the transmit powers based on the observed interference such that a node 
utility defined as a concave function of SINR minus a linear power cost, is maximized, i.e., distributed 
maximization of the Lagrangian for an optimization problem where, for arbitrary power costs (dual 
variables), the obtained transmit powers are not necessarily optimal. Other works, such as (H, investigate 
game theoretic approaches in power control converging to a Nash equilibrium. The analysis of convergence 
rates of some distributed power control techniques based on fixed-point iterations can be found in ifTOl 
and the references therein. 

A second category of research works investigate globally optimal power control in interference net- 
works, especially the works of Chiang et al. ifTTTl — lfT4l . The work in |[T4l covers some of the major 
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problems and the corresponding solutions. The research works in this category attempt to characterize 
the set of feasible transmit powers as a convex set and propose a convex optimization framework to 
find the optimal power levels for the wireless transmitters. One of the standard methods of describing 
the feasible set of transmit powers is to consider the spectral radius of a matrix made of network-wide 
parameters 0. Chiang et al. use this notion to propose a distributed power control algorithm in ifTTI - 
fT3l . Their work follows a network utility maximization (NUM) lfl5l framework where they use gradient 
ascend methods for a class of concave utility functions. While such algorithms are proven to converge 
to the globally optimal solution, the convergence speed is rather slow due to the reliance on the updates 
on the gradient of the Lagrangian. In addition, the algorithm involves finding the Perron root |[T6l of 
a matrix at each iteration. For an N x N matrix, finding the root is known to be of order 0(N 3 ) in 
complexity and therefore, not properly suited to larger network sizes. 

A different approach in solving network utility maximization problems is to formulate an equivalent 
problem whose feasible set is convex in its variables. For instance, it has been shown that the set of 
feasible power vectors is log-convex |[T4l . ifTTl , or, for a variety of utility functions, the optimization 
problem is convex if the problem is formulated in the inverse of power variables 1141 . Using the former 
technique, the work in lfl8l proposes a greedy power control technique based on the direction of the 
gradient of the objective function. As mentioned before, the major drawback of gradient-based methods 
is the potentially slow rate of convergence to the optimal solution. 

More recent works on power control focus on fast converging algorithms and non-convex cases not 
addressed before. For instance, the work in |[T9l proposes the MAPEL algorithm which solves the 
weighted sum-rate maximization problem. They express the objective function in terms of a product 
of exponentiated linear fraction functions; this leads to generalized linear fraction function programming 
(GLFP). However, there is, in general, a trade-off between the speed of convergence and global optimality 
in the solution |[T9l . Another work in EUl addresses the power control problem for a weighted sum-rate 
objective function. The authors propose iterative function evaluation for the solution and, empirically, 
show the convergence of their method. 

Regarding the speed of convergence in general resource allocation problems, recent works have pro- 
posed practical implementations of Newton's methods to attain the optimal solution with as little com- 
putational overhead as possible. Due to its reliance on fixed point iteration, Newton's method enjoys a 
quadratic convergence rate near the root. The work in ET1 investigates assignment of spectral resources 
in wireless cellular network where the authors improve the implementation of the Newton's method 
for the barrier function approach in solving convex optimization problems. In particular, for the system 
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model under consideration, the inverse of the Jacobian matrix can be calculated efficiently. Another work 
based on the same approach has been presented in ll22ll where, similar to EH, efficient calculation of the 
Newton's iterations leads to fast converging and numerically efficient resource allocators for OFDM-based 
cognitive radio networks. 

B. Our Contributions 

Motivated by recent works on fast resource allocation in cellular wireless networks, in this paper, we 
will develop a fast-converging power control algorithm. We start by formulating a utility maximization 
problem (containing NUM problems as a subclass) and argue that for a class of utility functions concave in 
the logarithm of SINR whose logarithm of gradient in terms of log-SINR is Lipschitz continuous |[23l . the 
problem can be solved efficiently (this includes most commonly used utility functions). We then propose 
an iterative algorithm based on fixed-point iterations, for problems with and without a minimum/maximum 
power constraint at the wireless transmitters. We prove that our proposed algorithm converges to the 
optimal power vector and the convergence rate is at least linear^. Moreover, for a network with N 
wireless transmitters, the complexity of our solution is 0(N 2 ). This is more efficient and scalable than 
gradient-based methods. In summary, our contributions are: 

• Formulation of the power control problem in an interference network for a general concave utility 
function. 

• Providing necessary and sufficient conditions for the optimality of the transmit power vector. 

• Presenting a fast and low-complexity algorithm which converges to the optimal solution for problems 
with or without power constraints. 

• Proof of convergence for the power control algorithm with or without power constraints. 

As will be illustrated in the numerical results section, one key benefit in using the proposed iterative power 
control is its capability in tracking slowly fading wireless channels. Moreover, since our formulation of 
the power control problem is based on a general utility function (as opposed to sum-utility), power control 
for more sophisticated network structures is possible, e.g., power control in cooperative cellular networks. 

C. Notation 

We use lowercase and uppercase bold letters to present column vectors and matrices, respectively. In 
line with this, (x)j (or simply Xj) represents the ith element of the column vector x and (A)jj (or simply 

'That is, at iteration k, the norm of error scales as 0(a k ) where < a < 1 
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Aij) is the element at row i and column j of matrix A. Furthermore, x T and A T are transposes of the 
vector x and matrix A, respectively. We use the function D(x) : 1^ — > M NxN to define the diagonal 
matrix with the elements of vector x as its diagonal elements, i.e., (D(x)) ii = (x)j. 

We extend the basic operations in field of real numbers R to the space 1^ as follows: For any two 
vectors x and y in WL N , we define the result of xy as a vector in whose ith element equals (x)j(y)j, 
i.e., the operation represents an element-wise multiplication of the two vectors. Similarly, we extend 
scalar functions of real numbers to the vectors in an element-wise fashion. For instance, (e x ) i = e^ x K 
and (x fe )i = (x)f for k G R. 

For any scalar function f(x) with x, f(x) G R, / G C 2 |^, f'(x) and f"(x) refer to the first and second 
derivatives at x. For a scalar function /(x) with xel" and /(x) G R, V x / G R^ and V x / G R NxN 
are the gradient vector and Hessian matrix of function / with respect to x. Subsequently, for a vector 
function f(x) : R^ -> R N , f'(x) G M NxN denotes its Jacobian matrix where (f'fc))^ = . 

Finally, for a square and symmetric matrix A, A >z and A ^ imples A is positive or negative 
semi-definite, respectively. For a given matrix A or vector x, A > and x > means all the elements 
of the matrix or vector are non-negative. 

D. Structure of the Paper 

The rest of this paper is organized as follow: Section Hl-Al presents our system model and the formulation 
of the optimal power control problem. It also presents the form of the optimal solution. In Section [Till 
we propose the power control algorithm achieving the power vector satisfying the optimality conditions 
presented in Section ITl-AI and provide sufficient conditions on the convergence of the proposed algorithm. 
Section [TV] is dedicated to numerical results based on a practical wireless network model, thereby 
validating the accuracy and efficiency of our proposed algorithm. Finally, Section [V] presents a summary 
of our contributions and concludes this paper. 

II. System Model 

Here we describe the system model under investigation. Then, we propose the utility maximization 
problem and the corresponding optimality conditions for the solution. 

2 C 2 indicates the class of functions with continuous second derivative. 
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A. Network Model 

We investigate a wireless network where N wireless contenders are communicating with their own 
destination over a shared medium, i.e., a logical network resource which is used by all wireless trans- 
mitters. We do not make any assumptions on uplink or downlink transmission or the method of channel 
access in the network, e.g., frequency /time duplexing or CDMA. Our methodology is based on the 
channels between each transmitter and all the receivers. We assume that vector h describes the power 
of the direct channel, i.e., hi is the channel magnitude squared for the link between wireless node i to 
its own destination. The interference channel is described by a matrix H where Hij is the (interference) 
channel magnitude squared from the link from node j to the destination of node i. Notice that we allow 
for self-interference, i.e., H^i > ^. 

To simplify the formulation, we define the normalized interference matrix V = D(h) _1 H. In this 
regard, the normalized received interference vector q (with qi as the normalized interference power 
received at the destination of node i) is defined as 

q = Vp + D(h)- 1 7 7 = Vp + C, (1) 

where p and rj are the transmit power and received noise power vectors (£ = D(h)~ 1 T/ represents the 
normalized noise term). Denoting the SINR vector at the receivers as 7, we have 

7 = -• (2) 

q 

Next, we will formulate the power control problem and obtain the corresponding optimality conditions. 

B. Problem Formulation 

Our goal is to find the optimal power control scheme maximizing a objective function of the SINR. 

4 



We consider an objective function U (7) : — > M to be a smooth, increasingly and concav43 function 

ogl U > 0, and Vi og ^ 



of the log-SINR vector , log 7 (that is, U G C 2 , Vi og7 f7 > 0, and Vw-y^ - °)- Based on tn i s objective 



3 As an example, phase noise at the receiver can be modeled as an additive term with a variance proportional to the received 
signal power. Otherwise, an assumption of Hi^ > leads to a finite cap on achieved SINR. 
4 We assume the utility of each link should increase with SINR. 
5 A concave utility function leads to fair assignment of network resources 1151 . 
6 We refer to such a function as log-concave. 
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function, the resource allocation problem can be formulated as 

max U(j) (3) 
7:P,q 

subject to: 7 = — (4) 

q 

q = Vp + C (5) 

Note that if U (7) = Yli=i ^t(7i) f° r concave increasing functions Ui(-), the optimization problem in ©- 
© turns into a network utility maximization (NUM) problem 031. The optimization problem in ©-([5]) 
is not convex due to the non-convexity of the constraint ©. However, we show that the problem can 
become convex in the logarithm of the optimization variables lfT3l and ifTTl . To this end, define 

x := log 7, y := log p, z := log q, (6) 

leading to the following equivalent optimization problem 

max U(e x ) (7) 
x,y,z 

subject to : z > log (Ve y + £) , (8) 

x < y - z. (9) 

It can be easily verified that the optimization problems in ©-(15]) and ©-© are equivalent(note that if 
y is a solution to ©-©, then the constraints © and © will be met with equality and therefore, © 
and © are satisfied). In order to prove the convexity of the optimization problem in ©, we require the 
following lemma: 

Lemma 1: The function /(x) = In (a T e x + a) where a S M. +N and a E IR + is convex in its domain 
(x G R N ). 

Proof: See Appendix lAl ■ 
We prove the convexity of the optimization problem © in the following theorem. 

Theorem 1: The optimization problem defined by objective function in © is convex optimization if 
and only if [/(•) satisfies the following 

D( 7 ) v' Un(i) + D( 7 )D(v T l7) ^ 0. (10) 

Proof: Note that from Lemma [T] the feasible set defined by © and © is convex. Now the problem 
is convex optimization if the utility function is concave in x, i.e., has a negative semi-defmite Hessian. 
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The condition for that is simply 

Vx^(e x ) = (Vxf/(e x ))' = (V)' T V 7 u)' (11) 
= D(e x ) y 7 f/D(e x ) + D(e x )D(v 7 £/) < 0. (12) 

where in (fTTT) the derivative (for the Jacobian) is with respect to x. By replacing e x with 7 in (|T2b . the 
requirement of the condition in ( fTOb is proven. ■ 

Note that the feasible set defined by constraints ©-© has none-empty interior and satisfies the 
qualification conditions for Slater's constraint |[24l and therefore, the duality gap between the optimization 
problem in ©-© and its corresponding dual problem is zero. Next we use the set of Kahrun-Kush- 
Tucker (KKT) optimality conditions to determine the optimal solution to the problem in ©-I©. The 
results are presented in the following theorem. 

Theorem 2: The optimal transmitter power vector p* maximizing the objective function in © satisfies 
the following system of nonlinear equations 

p* = P * v-rroCq*)- 1 (13) 

P P V T D( 7 *) y 7 t/Dtq*)- 1 V ; 

where 7* and q* are the corresponding SINR and interference power vectors for the given power vector 
P*- 

Proof: From KKT conditions, the gradient of the Lagrangian is zero at the optimal power vector. 
Following from ©-© we have 

L(x, y, z, A, n) =U (e x ) - A T (x - y + z) 

-/x T (ln(V e y + C)-z), (14) 

where A, /1 represent the Lagrange multipliers. From the gradient of L(-) with respect to x, y, and z we 
obtain 

y x i = =^=> D(e x )V 7 C7 = A (15) 
Vy L = 0^A T = /i r D( e y 1 VD(e y ) (16) 
Vz L = =► A r = fi T . (17) 

From ([B]) and CTJ) it follows for flU) that 

D(«-)V,I7 (18) 



V T D(e z ) _i D(e x ) y 7 U 
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where by replacing x, y, and z in (fT8T ) respectively with log 7, logp, and logq, the result in (fT3T > is 
obtained. We emphasize that the division in (fLSt is element-wise. ■ 
The optimal power vector is the solution to the system of N non-linear equations and N variables 
presented in Theorem [2] As explained in the introduction, we avoid using the Newton's method to 
solve (fT3T ) due to its computation cost for N S> 1. To this end, in Section UlTl we formulate an algorithm 
based on fixed-point iteration and prove its convergence for a wide class of concave utility functions 
[/(■). 

III. Power Control Algorithm 

Building upon the results of Theorem [2] in Section ITl-AI here we propose an efficient iterative algorithm 
to solve ([3]). We first consider the problem for the interference limited scenario, that is, we assume that 
compared to the interference level at the receivers, £ ~ and can be neglected from our formulatioi£ . 
Later on we will extend the results from the interference limited scenario to the scenario with £ > 
(which includes maximum and minimum transmit power constraints). 

A. Power Control under Interference Limited Circumstances 

Assuming £ = and defining S := D(7)V, from dH) and © it follows that 

Sp = p. (19) 

Note that p £ W N+ is the (right) eigenvector of the corresponding matrix S with an eigenvalue of 1. In 
the following lemma, we discuss the left and right eigenvectors of S and prove that p[S] = 1 where p[-\ 
represents the spectral radius of the matrix. 

Lemma 2: For any p > and S defined in (fT9l) , p is the unique eigenvector corresponding to the 
spectral radius of S where p(S) = 1. Furthermore, the left eigenvector of S, namely, a, corresponding 
to the spectral radius of S satisfies the following at the optimal power vector. 

a = D (q)- 1 V 7 U. (20) 
Proof: See Appendix iBl ■ 

7 Note that this assumption is equivalent to having no maximum power constraints on the wireless nodes. With no power 
constraint on the transmitters, the wireless nodes can transmit at very high power leading to interference levels far greater than 
the thermal noise power, i.e., relative to interference level we can assume £ ~ 0. 
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From Lemma [2] and (fT3l) we have the following equality at the optimal power vector 

P* = P* = PW), with 0(p) = (21) 

and a. defined according to (f2Qb (note that for an arbitrary p, a is not necessarily the left eigenvector of 
S). Equation (|2Tb suggests that the optimal power vectorp* is a fixed point of the vector function p0(p) 
(note that 4>(p*) = 1). This motivates us to propose a fixed-point iteration based on (|2TT) to solve the 
system of non-linear equations in Theorem [2] 

From Banach's theory on fixed-point iterations |[25l . the local convergence for smooth fixed-point 
functions relates to the spectral radius of the Jacobian matrix at the fixed-point. Specifically, if the 
Jacobian matrix has a spectral radius smaller than 1, the fixed-point iteration converge^. 

Before presenting the fixed-point iteration and proving the convergence of the algorithm, we require 
the Jacobian matrix of the vector function </>(•)• The following lemmas determine the Jacobian of </>(■) 
at the optimal power vector and some of its properties. To keep the formulation simple, for the rest of 
the paper, we define and use the vector r := £ as the ratio the power vector over a defined in (l20l) . 

Lemma 3: The Jacobian of vector function </>(•) at p = p* satisfies 

D(a)0'(p) = \s T (D(q)- 1 V* ^D(q)- 1 + 2D (r)" 1 ) S 

^(r^ + D^s)], (22) 

with: r = -, S := I - S, (23) 
a 

and I is as the N x N identity matrix. 

Proof: Ref. to Appendix ■ 
Lemma 4: The matrix 

D (r)» (S T D (r) -1 S - S T D (r)" 1 - D (r)' 1 S) D (r)3 , 

is negative semi-definite, has a spectral radius of 2, and has exactly one eigenvalue of corresponding 
to the eigenvector of (ap) 3 . 

Proof: See Appendix |Pl ■ 

8 In fact, Banach's theory expresses the convergence conditions in terms of the norm of the Jacobian. However, for a matrix 
M and any e > 0, there is a definition of norm || • || such that p(M) <|| M ||< p(M) + e (See 1261 ). Therefore, the condition 
can be extended to the spectral radius instead. 
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Lemma 5: Assume the log-concave utility function U(-) satisfies the following^ 

7 



D 



Then, at the optimal power vector p = p*, we have 



< B, V7 > 0. 



(24) 



D(r)= D(a)0 / (p)D(r)t ^ 0, 
D (r)5 D (a) 0'(p)D (r)*] = p [D(p)0'(p)] =< AB - 2. 

Proof: See Appendix El ■ 
Now we can propose the following fixed point iterations which converges to the optimal solution of ([3]) 
for a class of utility functions described in Lemma [5] 

Theorem 3: Consider the optimization problem © with a smooth log-concave utility function U (7) 
satisfying d24b . Then, V 6 € (0, Q max ), with # maa; = 2B 1 _ 1 , the following iterations converges to (an) 
optimal solution of ([3]) 



Pk = 4>e(Pk-i) = 0Pk-i<f>(Pk-i) + (! - 6>)p fe _i, 

with arbitrary po G 



(25) 



Proof: The proof is complete if we show that at the optimal point p, p [0g(p)] < 1 ll25l . |[26l . It 
simply follows that 



p[(0 e (p))'] =p[eD(p)0'(p) + l] 



6»D (r)5 D (a) 0'(p)D (r)= + I 



(26) 



where d26l) readily follows from the definition of r in Lemma From Lemma [5] we know that the 
matrix in (l26l ) has a largest positive eigenvalue of 1 and a largest (in magnitude) negative eigenvalue of 
at most 4(B — 1). Therefore any 9 < 2 b~\ guarantees that p[<j>g(jp)'] < 1. On the other hand, from 
Lemmas |4] and |5] we know that for matrix <p'(p), p is the only eigenvector for the only eigenvalue of 1 
, i.e., p[4>' g (p)] = 1. However, this does not contradict with the convergence conditions since p is the 
optimal power vector. The reason is that if at iteration k, Pk is close to the optimal power vector, then 
the projection of the error onto the direction of the optimal power vector remains intact while the other 



9 Equivalently, log y 7 f/ is Lipschitz continuous in log 7 with a Lipschitz constant B. 

10 Note that for any two matrices A and B where B is invertible, p[A] — p[B _1 AB] and p[AB] = p[BA] 
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components of the error are attenuated based on the other eigenvalues of (/>(p)rH Hence, the algorithm 
converges to the optimal power vector p (or a positive multiple of p). ■ 

From theorem (0) we realize that while the iterations converge to the optimal power vector, the limit 
is not unique. In the following corollary we propose a modified fixed-point iteration which converges to 
a unique optimal power vector with finite norm. 

Corollary 1: The following fixed point iteration converges to a solution p for the problem in ©-(O 
where || p ||2= 1. 

Pfc = <t>6 ( j. 1 • (27) 



Pfc-l ||2, 

Proof: See Appendix 10 ■ 
The algorithm proposed in Theorem [3] requires iterative evaluations of function </>(•)> which requires 
0{N 2 ) calculations. Therefore, the fixed-point iterations based on (f2Tb will have a complexity of 0(N 2 ). 
We would like to emphasize that this level of complexity is less than both gradient-based and Newton's 
iterations proposed in the literature and as a result, leads to faster convergence for N S> 1. In Section lTlI-BI 
we extend the results of Theorem [3] to the scenario with non-zero additive noise and show that the 
algorithm converges for an even larger set of utility functions. 

B. Power Control with non-zero Noise Power 

In this scenario, we assume £ > 0, and consider minimum and maximum power constraints at the 
wireless transmitters, i.e., for wireless user i, p m in,i < pi < p m ax,i (Remember assuming non-zero noise 
power but no power constraint at the transmitters leads to the interference limited scenario). We start by 
deriving the optimality conditions for a power constrained problem. We prove the results for the case 
where each wireless transmitter has a maximum transmit power constraint. The analysis for the scenario 
with minimum power constraints is very similar. 

In order to proceed, we need to add a new constraint to the optimization problem in {/J)-©, i.e., 
y = lnp < lnp max . In the following theorem, we provide the necessary and sufficient conditions for 
the optimality of a solution to the power-constrained problem. 

Theorem 4: The optimal power vector p* , maximizing the optimization problem Q under the transmit 
power constraint p < p ma x, satisfies the following equation 

p* = min {pmax, P*<f>(p*)) • (28) 
"Also, note that in the interference above thermal scenario, if p is an optimal solution, then ap for a > is also optimal, 



.e., the optimal solution is a line in R + and is not unique. 
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Proof: Without loss of generality assume p m ax,i = 1 for all nodes (for p m ax,i ^ 1> normalize row 
j of V and also Q to p m ax,i), leading to the constraint y < 0. Since the new constraint is convex, we 
can use the KKT conditions to identify the optimal solution. Define the new Lagrangian to be L c (-) and 
consider [3 > as the the new dual variable. The gradient of Lagrangian with respect to x and z is 
similar to ( fT5T ) and ( TTVT ) respectively. From the gradient of L c {-) with respect to y it follows that 

Vy L c = => X T - [3 T = A T D (e 2 )- 1 VD (e x ) (29) 

where after replacing x, y, and z respectively with 7, p, and q we obtain 

P* = (30) 

From above, if = then p* = (<p(p*)i) p*. If > 0, by the optimality conditions we must have 
y* = which leads to p* = 1. On the other hand, a* -f3*/p* < a* and therefore, p*(0(p*))j > 1 which 
satisfies p* = min(l,p* (^(p*))^. ■ 
For the problem with minimum constraints on the transmit power of wireless nodes a similar result can 
be obtained. 

Theorem 5: The solution to the optimization power control problem presented in © with the maximum 
and minimum power constraint (p m j n < p < p m ax) satisfies the following 

p* = max(p m j n ,min(p mo:r ,p*0(p*))) . (31) 

Proof: The sketch of the proof is similar to that of Theorem |4] and can be easily repeated by the 
reader. ■ 
The power control algorithm converging to the optimal solution of (f3]) under maximum power con- 
straints is similar to that of interference limited networks. Before proceeding with the proof of conver- 
gence, we present the Jacobian of 4>e(-) (defined in (T25T )) at the optimal power vector in the following 
lemma. 

Lemma 6: For the optimization problem in © with £ > and maximum power constraints, the 
Jacobian matrix of 4>e{-) at the optimal solution satisfies 

0' fl (p) = D(0(p))0£ >c=o (p) + J, (32) 

where cf>' g ^_ (p) is the Jacobian at the optimal solution as if C = and J is a matrix with the property 
that for any vEl N and any i that (0(p)) i = 1, (Jv) i = 0. 

Proof: See Appendix iGl ■ 
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In the next theorem, we present the fixed-point iteration that converges to the optimal power vector 
under maximum transmit power constraints at the wireless transmitters. 

Theorem 6: The following fixed-point iteration converges to the global optimum of power control 
problem defined in Theorem [3] under maximum power constraints at each node, 

p fc = min (pmax, <j>e(Pk-l)) = mm (1, e (p fc _i)) . (33) 
Proof: First note that from Theorem HI any fixed point of (l33l is an optimal solution to the 
optimization problem defined in © under maximum transmit power constraints. Since the function in (l33l 
is not differentiable at the optimal power vector, the proof is slightly different from that of Theorem [3] 
Without loss of generality assume {4>e{p))i > 1 ( or equivalently, (0(p))j > 1 ) for i € {1,2, ...,n} 
and (4>e{p))i = 1 for i £ {ra + 1, n + 2, ...,7V}. We know that (f>(p) is continuous in its domain and 
therefore, for any e > 0, we can find r > small enough such that sup p / gB ( p r ) || 4>e(p') ~ 4>e(p) ||< e 
where B(p,r) is an open ball of radius r (in R^) centered at p. Choose the r corresponding to an 
e < mini< i < ri (</> e (p)) i . 

Now suppose at iteration k we have p^ = p + £fc with e k € 73(0, r). Define A„ as a diagonal matrix 
with (A n )j j = for i < n and (A n )j j = 1 for i > n. It follows for the power vector at iteration k + 1 
that 

Pfc+i = min (1, e (p + e k )) ~ (I - A n ) 1 

+ A n (p + e/>e(p)e fc ) = p + A n <p' 9 {p)e k . (34) 

Note that (l34b comes from the continuity and linear approximation of </>e( - ) in the vicinity of p. From (l34l) 
and Lemma [6l we have 

efc+i = Pfc+i - P - A n 0g(p)e fc = A n </) / 0)C=o (p)€ fc , (35) 
where (|35T ) is true since A n J = OatxA^ and A n D (</)(p)) = A n . Furthermore, from the proof of Lemma[6] 
we have p A n cf)' e ^ =0 (p) < p 4>g ^=o(p) < ^ an( ^ therefore || e^+x ||< p || || for some consistent 
norm function || • ||. This concludes the proof of local convergence of the iterations. ■ 
Based on a similar approach, we can propose the following fixed point iterations for the scenario with 
both maximum and minimum power constraints. 

Theorem 7: For any U(-) and 6 satisfying the conditions of Theorem [3l the following iteration con- 
verges to the optimal solution of problem ([3]) with the minimum and maximum power constraints p m in 

and Pmax 

p k = max (p 

min-, 

max (p 

max 5 ^e(Pfe-i))) • 
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Proof: The proof is similar to that of Theorem [6] and is not repeated here due to space constraints. 

■ 

The convergence of fixed-point iterations under maximum and minimum power constraints at the wireless 
transmitters depends on the stability of the Jacobian matrix 0# £_ (p) which is, in turn, guaranteed if 
the utility function satisfies the conditions mentioned in Theorem [3] and 9 is chosen small enough. 
The following corollary extends the proof of convergence to all log-concave utility functions for an 
optimization of the form of ([3])- © under minimum and maximum power constraint. 

Corollary 2: For any strictly increasing log-concave utility function U(-), the fixed point iteration 
defined in Theorem [7] converges to a point which satisfies the KKT conditions for the optimization 
problem © under maximum and minimum transmit power constraints. 

Proof: Since the utility function is smooth, over compact set S = {p : p m in < P < Pmai}. we have 
05/ (p) < Bg for some finite positive B$. Then, for any 6 < 8 max < 2 bI-i t ' ie P ower control algorithm 
is locally convergent. ■ 

C. Asynchronous Power Updates 

As a final note, a property that may be of implementation benefits is that the proposed power control 
algorithm does not require synchronous power update^. In fact, the convergence speed might be faster 
under asynchronous power updates. The proof of convergence under asynchronous power updates follows 
a general property of fixed point iterations. Consider x = f(x) a converging iteration to a point xo with 
p [f'(xo)] < 1. The asynchronous power update process can be formulated as follows 

Xfc+l = gJv(gJV-l(-(gl(Xfc))...)) 

with: 

gi (x)=D(e i )f(x) + (I-D(e i ))x (36) 

and ej is a vector with a single non-zero element (e,) i = 1. The new fixed point iteration defined in (|36l ) 
converges if the Jacobian matrix is bounded by 1. It follows that 

(gjv(gJV-l(-(gl(x )). ..)))' = g7v( x o)gjV-l(x )-gi(x ) 

N 

= H (D( ei )f'(x ) + I - D(eO) = f'(x ), 

i=l 

suggesting that the asynchronous methods also converges. 

^Asynchronous updates is defined as the process of calculating the new value of variable (xk+i)i based on (xfc)^ for j > i, 
and (Xfc+i)^ for j > i. 
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D. Distributed Power Control and Implementation Issues 

In order to solve the optimization problem ([7]), the destination nodes need to know the interference 
matrix. This represents a significant signaling overhead of 0{N 2 ) which might not be practical in dense 
wireless networks. Here we show that through distributed implementation, the signaling overhead can 
be reduced to 0{N). From the fixed-point iterations in Theorems [3] and [6] it is observed that at a given 
iteration, the destination of wireless transmitter i requires the values of aj and (S T a). to find the new 
transmit power. On the other hand from Lemma |2j the value of aj depends on the gradient of the 
utility function as well as the received interference power, q t . Assuming the receiver can estimate the 
value of qi, knowledge of {sj-yU)i is required for the evaluation of Oj. In order to evaluate (V7^(7))i> 
global knowledge of the SINR vector 7 is necessary. Therefore, we assume that at the beginning of each 
iteration, every destination node broadcasts the SINR for the corresponding wireless transmitter. 

Furthermore, for (S t q) it follows that 

N N N 

(s 7 '*), = Y s '-J n J = Y v J- r '°j = Y 7/ <•''•> 
j=i j=i j=i 

where K{ = Therefore, in order to evaluate (S T a) . at each destination node, global knowledge 
of vector n at each iteration, and local knowledge of channel gains is required. Note that this local 
knowledge of channel gains includes channels gains from wireless transmitter to other destination nodes 
in the network which can be provided through channel feedback or channel reciprocity assumption in 
the case of TDD channel access. 

Algorithm 1 provides a distributed implementation of the proposed power control and converges to 
the optimal solution of ©. Note that implementation of Algorithm 1 requires 2N messages broadcasted 
at each iteration. More importantly, evaluation of the power vector depends on the local channel gains 
which leads to a signaling complexity of O(N) in between the destination nodes. From Algorithm 1, 
two phases of signaling are required at each iteration: During the first round the destination node(s) 
broadcast the SINR vector to allow other destinations to evaluate the gradient of the utility function, and 
subsequently, the value of ex. The second round involves broadcasting the value of ex. This information 
is used to determine how much to increase/decrease the transmit power for the next iteration. Note that 
an assumption of NUM for the optimization problem obviates the need for the first round of signaling 
(each node can find the derivative of its own utility individually based on its own SINR and interference 
power). This suggests an overall N signaling messages per iteration. Finally, in Algorithm 1, the value 
of 6 is adaptive, that is, after each M iterations, 8 becomes smaller. In general, the initial value of 9 and 
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1: Variables at user i, iteration k: p\ k \ q~ k \ 7 t - , af 1 , <f>\ , and k\ . 
2: Initialize pf ] G (0,p ma x], G (0, 1], M G N, and k,c = 1 
3: repeat 

(k) 

4: Estimate received interference power q\ 

(k) n( fe ) 

5: Broadcast j> ' =«- 

nil, ( fc ) , (v*V( 7 «)) 

6: Calculate a) ' <— - m — Ll 

1 q. ' 

7: Broadcast nf> <- 7 ' 
8: Calculate </> (fc) ' ^ 



9: Update power: p\ <- mm (p ma x,P- <p\ 
10: Calculate error : e ( fe+1 ) <-| pf +1) -pj fc) | 
ll: if A; > cM then 

12: <- § 

13: C^-C+l 

14: end if 

15: k k + 1 
16: until e» +1 ) < tol. 

Algorithm 1 - Distributed Power Control: General Utility Function 
M could be tweaked to obtain better convergence properties than that of fixed 9. 

IV. Numerical Results 

In order to validate the theory developed and illustrate the performance of our proposed power control 
algorithm we present the results of simulations based on popular system models. We consider a seven-cell 
cellular topology with hexagonal cells where, within each cell, N ce u wireless users are connected to the 
central base station. We consider a dense interference network where, due to the short distance between 
the base-stations, inter-cell interference imposes significant performance degradation. We adopted the 
IEEE 802.16-based channel models for fixed wireless network with intermediate terrain and path loss 
conditions ll27l . The simulation parameters are presented in Table J] Fig. Q] illustrates the distribution of 
wireless users across a 7-cell network. 

In the examples, we consider the log-rate as the utility function, i.e., 17(7) = logi?( 7 ) where R(-) is 
the achievable rate for a SINR of 7. This utility function leads to an efficient balance between overall 
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Parameter 


Value 


1>I UlllUCl Ul LClln 


7 


1>I UlllUCl Ul Uacla pel UC11 




Channel Access 




\^.aillCl ncLjuciiL-y 


±.U V.J 1 1 /. 


\^.Cll ivalllUS 


ouum 


A trfpntr^i IVToHpI 

/T-llLClllltl 1V1UU.C1 


OMNT - 15 HR 


Path Loss Exponent 


3.79 


Shadowing 


Log-normal STD = 9 dB 


Channel Coherence Time 


10ms, oo 


Power Update Intervals 


5ms 



TABLE I: Simulation parameters. 
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Fig. 1 : Distribution of wireless node in a 7-cell network topology. 

throughput and fairness according to Jain's fairness index (weighted proportional fairness). Assuming 
R = log 2 (l + 7), we have —1 < -^jtt^ < 0. Then, we have B = 2 (see Theorem [7]) and any 
6 < &max = 1.0 will guarantee convergence. 

A. Rate of Convergence 

In our first example, we consider a 7-cell static wireless network with parameters chosen according to 
Table U We assume random channel instances and initialize the power control problem with a random 
starting point every time. We set R(j) = log(l + ^) corresponding to a 7 dB gap to capacity. Fig. [2] 
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5 10 15 20 25 30 

Number of Iterations 



Fig. 2: CDF of number iterations required to converge to 5% and 1% of optimal solution. 

plots the cumulative distribution function (CDF) of the number of iterations required to get to 5% and 
2% vicinity of the optimal solution. Based on the results, in 90% of instances, the algorithm converges 
within 10-15 iterations. 

B. Robustness to Channel Variations 

Thanks to the fixed-point iteration nature of our proposed algorithm, the convergence speed is at least 
linear. Therefore, in general, the algorithm can track small changes to the system parameters in a very 
short time. To this end, we consider the same topology as in the previous example, except that this time 
we include slow channel fluctuations. Specifically, we assume that the channels are constant within a 
coherence time of 10 ms and change independently. 

We consider three scenarios: first, the optimal case where the channel powers are updated instanta- 
neously; second where the power update period is half the channel coherence time, and third, where the 
only channel information available is the average path loss. The results of the simulation are presented 
in Fig. [3] As is clear from the figure, the algorithm can efficiently track the channel variations. 

C. General Utility Functions 

In the previous two examples, each node had its own utility function, i.e., there were illustrations of 
network utility maximization . However, our theoretical development is valid for a broader class of utility 
functions. In this example we illustrate the ability of the algorithm to optimize power in networks where 
the utility function is coupled across nodes. We consider a simulation scenario with relay nodes added to 
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Fig. 3: Robustness of the power control with respect to small fading channel correlation in between power 
updates. 

our cellular network. We assume in each cell, there is a dedicated relay which helps wireless transmitters 
with poor channels to the base station, e.g., cell-edge users. The relay operates in decode-and-forward 
mode ll28l . i.e., the transmission time is divided into two equal time slots where during the first slot 
the wireless user sends its message to the relay, and during the second time slot, the relay decodes and 
re-encodes the user's message and forwards forwards it to the base station. It follows that the achievable 
rate under cooperation is given by 

1 / min(7 w ,7 r6 )\ 
R c = - log 1 + , (37) 



2 ° v r 

where ^ wr and ^ are the corresponding SINRs for the wireless transmitter-relay and relay-base station 
links and T > 1 is the decoding loss. The function min(-, •) has a discontinuous derivative and therefore, 
we use the following analytic approximation for the power control algorithm 

mm(x, y) ~ ^— log (e~ kx + e~ ky ^j , k > 1. (38) 

Choosing a lager value of k in (I38T ) leads to better approximation of min(-,-) but at the same time 
increases B defined Theorem[3 leading to slower convergence speed (in our simulations, we select k = 5 
and 6 = 0.5 for a log utility 1 13 L Note that the min(-, •) introduces a coupling of power across the relay 



and source nodes and the Hessian of the utility function is not diagonal. This is due to the achievable rate 
13 Furthermore, the logarithm of rate based on the proposed approximation is concave over the range x,y > t. 
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for the wireless users using a relay, which is a function two SINRs, namely, user-relay and relay-base 
station SINRs. 

For a wireless transmitter not using the relay, the achievable rate is 



where 7^.1 and j s d,2 are the corresponding SINRs in the first and second times slots (the sources of 
interference are different in the two times slots if relays are used in adjacent cells - the source in the first 
time slot and the relay in the second). Unfortunately, using d39l ) in the optimization problem © leads 
to a non-convex problem. Therefore, we use the following upper-bound on the achievable rate for the 
optimization problem 



In other words, the upper bound in d40l is the capacity of the link to the base station when the interference 
is constant and equal to the average interference received from other cells (from wireless users and relays). 
Fig. |4] illustrates the achieved sum utility in this network for different power control and transmission 
techniques. In the case of power control, the algorithm is initialized with maximum transmit power and 
when no power control is used, the wireless transmitters transmit at maximum power. Based on the 
results in Fig. 01 relaying helps with the overall quality of service in cellular networks. Moreover, there 
is considerable gain in power control after relaying which made feasible through applying the power 
control algorithm presented in this work. Note that the difference between the convergence speeds of the 
relayed and non-relayed scenarios is due to the difference in the value of 9. Smaller 9 assures stability 
at the cost of convergence speed (however, for any 9, the local convergence rate is still linear). 



In this paper, we considered the power control problem in a dense interference environment for a 
general utility function. We proposed the form of the optimal power vector as the roots of the KKT 
equations. By introducing an algorithm based on fixed-point iteration, we solved for the power vector 
and showed that the proposed algorithm converges for a variety of log-concave utility functions. We 
considered both power constrained and unconstrained scenarios. In particular, we proved that for the 
scenario with both maximum and minimum power constraints on the wireless transmitters, the proposed 
algorithm attains the optimal solution for any given log-concave utility function, including functions that 
do not lead to network utility maximization. 




(39) 




(40) 



V. Conclusion 
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Fig. 4: Power control in a relay-assisted cellular network. 

Our numerical results shows that the algorithm converges in a few iterations. Furthermore, due to 
the fast convergence, the algorithm is very robust to channel variations. Finally, it was shown that the 
proposed algorithm is capable of solving a general class of utility maximization problems. As an example, 
power control in a relay-assisted cellular network was demonstrated. 

Appendix A 
Proof of Lemma Q] 

The proof shows that the Hessian of /(•) is positive-definite. It readily follows from the definition of 
/(•) that 

Vx/ = (Vx/)'=(^D(a)e^ = - ?; l T -x 

f 1 " „ri ■ D( a ) D ( fiX ) llT ) D ( a ) D ( eX ) • < 41 ) 
Define u = D(a)e x , i.e., a T e x = l T u. Then it follows for the Hessian of /(•) that 

(l T u + a) 2 Vx/= (l T u + a)D(u) - uu T . (42) 

Since u > 0, we can multiplying the sides of (l42l by have D (u) 5 y- (Note that this conserves the 
positive semi-defmiteness). Then, after multiplication, the right hand side of (l42l) equals 

(l T u + a)I-D(u)hl T D (u)l (43) 
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Note that the second matrix on the right hand side of (|43T > has a single non-zero eigenvalue of || 

i „ 

D(u) 2 1 ||2> and the rest of the eigenvalues are all zero. On the other hand, it simply follows 
that || D (11)2 1 || 2= || u || 1= l r u. This suggests that the eigenvalues of the matrix on the right hand 
side of (03]) are all non-negative and therefore, /(•) is convex. ■ 

Appendix B 
Proof of Lemma |2] 

Here we present a short version of the proof. For a more comprehensive study of such matrices we refer 
the reader to 0, lfT2l or any text on Perron-Forbenius norms of a positive matrix. Note that for matrix S 
we have Sij > Vi, j. From Perron-Forbenius theory, p(S) corresponds to a positive eigenvalue, namely, 
A. It suffices to show that A = 1. Note that if A > 1 then there is a vector p' > where Sp' = Ap'. 
This suggests that D(7)q' = Ap' where q' = Vp' is the interference vector corresponding to the power 
vector p'. On the other hand we know that p' is the eigenvector corresponding to the eigenvalue of 1 
for the matrix D(7')V where 7' = E7 which implies 7' = A7. Therefore, 

D(V)Vp' = AD(7)Vp = A 2 p'. (44) 

This result contradicts with the fact that p' is the eigenvector corresponding to the eigenvalue of 1 of the 
matrix D(7')V. Therefore p[S] = 1. For the second part of the lemma it simply follows from (fT3l that 

= D( 7 ) y 7 U = D(q)- 1 y 7 y 
V^D(7)v 7 ^D(q)- 1 P V^D(7)v 7 ^D(q)- 1 
where dividing the sides by the non-zero vector p yields, 

V T D( 7 ) y 7 ^D(q)- 1 =D(q)- 1 V 7 C7, (46) 

i.e., x/^UYi (q) 1 is an eigenvector of S T = V T D(7) corresponding to its largest eigenvalue of 1. By 
the Perron-Forbenius theorem, this eigenvector is unique and, therefore, it has to be a. ■ 

Appendix C 
Proof of Lemma [3] 

Consider p and a to be the corresponding right and left eigenvectors of S at the optimal point. Note 
that at the optimal power vector p, we have (f>(p) = 1. Then it follows from (f2TT > 

0/(p) = (s^) , = D (( sTa )") Q/ 

-D(a(S T c t )" 2 ) (S r «)' 

= D(a)- 1 (a' - (S T a)') . (47) 
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Following (S T a.y we have 



(S T a)' = (V T D( 7 )a)' = V T (D( 7 )a' + D(a)V) 
= S T a + V T D(a) 7 '. 



For 7' it follows that 



Y=(E)WV-D(£) q ' 

= D(q)- 1 (I - D( 7 )V) = D(q)- 1 (I - S) , 
where we use the fact that q' = V. Finally, for a' and from Lemma |2] we have 

o , = (D(q)- 1 V 1 CO' 

= D(q)- 1 (v 7 f/)V-D^) V 
= D(q)- 1 V * ^D(q)- 1 (I - S) - D (™) D( 7 )V 
= D(q)~ 1 V * t/D^S-D^S. 



Note that in derivations above, S = I — S, a = . Assuming r = we can use (|48l ), (|49l ), (l50T > 



(48) 



(49) 



(50) 



in (|47T ) to obtain 



0'(p) = D(a)" 1 [^(q)" 1 V ' ^(q)- 1 ^ 
-S T D (r)" 1 S - V T D (r) _1 S 



D(a)" 1 ^(q)- 1 V ^D(q)- 1 S 



S T D (r) _1 S - S r D (r) _1 S 



1-1 



(51) 



Finally, by inserting I— S in (l5Tb instead of S, and multiplying the sides from left by D(a), the expression 
in Lemma ([3]) is achieved. ■ 
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Appendix D 
Proof of Lemma @] 
It follows for the expression in the lemma that that 

D (r)5 S T D (r)^" D (r)^ SD (r)s 
- D (r)3 S T D (r)^ - D (r)^ SD (r)3 
= (D (r)5 S T D (r)^" - i) (d (r)^ 1 SD (r)» - I! 

D (r)» S T D (r)^) (d (r)^ SD (r)*) - I. (52) 



Note that all the elements of the first term in (1521 ) are positive. From Perron-Forbenius theorem, there 
is a unique positive eigenvector corresponding to the largest eigenvalue (spectral radius in this case). 
By observation, we realize that this eigenvector is simply (ap) 2 which leads to an eigenvalue of 1. By 
uniqueness, we realize that (l52l is negative semi-definite with exactly one eigenvalue of corresponding 
to the eigenvector (ap) 2 . ■ 

Appendix E 
Proof of Lemma [5] 

It follows for the expression in the lemma that 

D (r)» D(a)0'(p)D (r)» = D (r)> [s T (D(q)- 1 V ^ U 
D(q)" 1 + D(r) _1 ) S D(r) 2 +D(r) 2 S T D (r) -1 S 

-S T D(r) _1 -D(r) _1 sl D(r)* . (53) 

First note that by Lemma HI the second term in (l53l is negative semi-definite with a spectral radius of 
2. On the other hand, the first term depends on the Hessian of the utility function with respect to 7. It 
follows for the Hessian that 

DW - Mrow -> = D (f)* D (E)* D (.L) i 

^ ro (D' D (=0* D (f)' 

=D(r) " D (^) Iv ^ D (^f) ID(r) " 

= D(r)^M c/ D(r)^, (54) 
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where M[/ := D 
we have 



^-£7 i V7 UY) ( ) * . On the other hand, from (HOI) and the assumption in (1241 ). 



-(5-1)1 ^Mjj + I^O. 
By using (l55T ) in the first terms of (I53T ) it follows that 
D fr)5 



(55) 



DM 



S r (D(q)" 1 V * fZD(q)- 1 + D (r)" 1 J S 
= D (r)» S T D(r)^ (M v + I) D(r)^SD (r)» (56) 
b -(B-l)D(r)»S r D(r) -1 SD(r)5 . (57) 
From (|55T ) and (|56"1) . we see that the first term in (l53l) is negative semi-definite. Furthermore, following (|57T ). 



we realize that 



D(r)sS J D(r) _1 SD (r) 



where (1581 ) comes from the properties of matrix normo It follows that 



p(S T S 



gl.2 



2- 



(58) 



S lU = max II Sx 1 1 2 = max II x — Sx IU 

||x|| 2 =l l|x|| 2 =l 

< max || x H2 + || Sx IU 

||x|| 2 =l 

< max || x H2 +p(S) || x ||2< 2, 

ll x l|2 = l 



(59) 



where the last result is an strict inequality since x cannot be p (in such a case, || Sp ||= 0). Finally, 
from (I53T ). (l57l ) and (l59l we realize that D (r)5 D(a)0'(p)D (r) 3 is negative semi-definite and also 



p(D(p)0'(p)) < 4(5-1) + 2. 



(60) 



Appendix F 
Proof of Corollary [I] 

First note that p is a fixed point of d2"7T ). More importantly, for any r > 0, rp is «of a fixed point 
of (127T ). From the Jacobian of (1271) it follows that 

p \Y ±1 1 



4>o 



p II2 



p II2 



p II2 



(I + 0D (r) A) (I- PP T ), 



(61) 



14 For any positive semi-definite matrix A, we have p(A T A) 
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where A = D(a)0'(p) is a symmetric negative semi-definite matrix with a single eigenvalue of 
corresponding to the eigenvector p (See Lemma 0. From linear algebra we know that the eigenvectors 
of J are linearly independent and form a basis in R^. Consider Vi to v^r as the eigenvectors of A 
corresponding to the eigenvalues Ai to A at where Ai = 0, and A^ < Aj_i are all real and negative. 
Furthermore, A at > — 4B + 2 where B is defined in (O. 

We need to show that the spectral radius of (|6"TT) is less than 1. It follows that 



P 2 



p[(l + ffD (r)A) (I- PP T )] 

'i + a)q , 



where 



A = D (r)* AD (r)~ 
Q = I - D (r)= pp T D (r)^ 



(62) 

(63) 
(64) 



Note that A and A share the same eigenvalues. Furthermore, D (r x ) 2 Vj is the eigenvectors of J 
coiTesponding to the eigenvalue Ai. By assumption, we have vi = p and since QD (r-- 1 ) 2 p = 0. For 
the other eigenvectors Vj, 1 < i < N it yields 



1 + A) QD (r)~ v, = [1 + Aj D(r)~ Vj 

= (1 + 0Aj)D (r)^ Vj, 



(65) 



which readily suggests 



P 



P 2 



max I 1 + 9Xi |< 1.1 

,2<i<N 



Appendix G 
Proof of Lemma [6] 



First note that for 4>'(p) we have 

0'(p) = D 
>(P 



q'-D 



a 



a. 



D 



) (a'-(I + D(0(p))-I) (S T «)' 
D(0(p)) U =0 (P) ~ (D(0( P )) - I) (S T a)' 



(66) 
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where <^ =0 (p) is the Jacobian of <fi(p) in the interference limited case evaluated at p. For <fie{p)' it 
follows that 



Note that (fie £=o(p) = I + #D(p)0^ =o (p) is simply the Jacobian of <fie(p) assuming C = 0- The only 
difference is that, since p is no more an eigenvector of </>6»,c=o(p) corresponding to an eigenvalue of 1, 
p \4>e £=o(p)] is strictly less than 1. The second terms in ( f67T ) is also a bounded matrix. Also note that 
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